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Abstract 

Phase transition of decay rate from quantum tunneling to thermal activ- 
ity regimes is investigated in (3+l)-dimensional field theories with symmetry 
breaking term fcft. By applying the two independent criteria for the sharp 
first-order transition to the same model, the upper and lower bounds of crit- 
ical value of the symmetry breaking parameter are obtained. Unlike two 
dimensional case continuum states of the fluctuation operator near sphaleron 
solution play an important role to determine the type of transition. 



Recently, much attention is paid for the phase transition of the decay rate from thermal 
activity to quantum tunneling regimes in (l+d)-dimensional field theories when d = 1[1], 
d = 2 [2], and d = 3 [3]. Especially, d = 3 case, which is a concern of the present letter, has an 
abundant applications in cosmology[4] and particle physics[5, 6]. In this problem potential 
is usually choosed, for simplicity, as an asymmetric double well potential 

W] = ^(0 2 -i) 2 -M (i) 

where / is symmetry breaking parameter whose range is < / < 4/3>/3. In fact, Eq.(l) 
is a tree level potential of the temperature-dependent corrections[7,8]. The analysis of the 
phase transition for the temperature-dependent effective potential by the method used in 
this letter is expected to be very complicated and will be discussed elsewhere. After taking 
an appropriate scaling transformation it is easy to show that the system with potential (1) 
exhibits a same physics with a system whose potential is [9] 

VM = \<P A -<P*+ l -5<P\ (2) 



if 5 has a relation with / as follows: 



/=^4. (3) 
(3-6)* 
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From the range of / one can see that the parameter 5 can vary between and 2. In the limit 
5 — > 2(or / — > 0), we can use the usual thin wall approximation method[10]. In the opposite 
limit 5 — > 0(or / — > 4/3 a/3), the height of the barrier between the two vacua vanishes. Since 
the radius of bubble is smaller than the wall thickness, this limit is frequently called thick 
wall limit [4]. 

In Ref.[3] Ferrera has shown that the type of the phase transition of decay rate is sharp 
first-order when / = 0.25 and / = 0.55, and smooth second-order when / = 0.75 by 
solving the equation of motion numerically with an aid of the multigrid method. This means 
that there exists a critical value /*(or <5*), which distinguishes the type of transition. The 
multigrid method, however, is not useful to determine /*(or 5*) because it is impossible to 
treat temperature continuously due to the fixed grids. Furthermore, it may be highly difficult 
to solve the equation of motion numerically in realistic models. It is, therefore, important 
to determine /* without solving equation of motion. It may be worthwhile noting that the 



terminology "phase transition" used in the present letter does not mean the generic transition 
of the physical system's phase but the transition between the two different instanton regimes. 

In this letter we will evaluate the lower and upper bounds of /*(or 5*) by using two 
independent criteria for the sharp transition. 

Let us start with Euclidean action 
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which yields an equation of motion 



d 2 



dr 2 



Since zero temperature bounce solution (fts has 0(4) symmetry[ll], the equation of motion 
for <pB is reduced to 

d 2 (j)B 3 d(f>B 



+ 



V' 



2l9B\ 



dR 2 R dR 

where R = \/t 2 + x 2 . Thus, the classical action for 0(4) symmetric solution is 



(5) 



(6) 



At high temperature thermal activity is governed by static sphaleron solution [12] which 
has a 0(3) symmetry. Hence, the equation of motion for the sphaleron solution <f> sp h is 

d 2 4> sp h 2 d(j) sp h 



+ -- 



dr 2 r dr 

where r = y/ii?, and, thus, the classical action for this configuration is 

S S ph — TsphSs- 

Here, r sp h is the period of small fluctuation around (fi sp h and 

2 



(7) 



(8) 
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= 4tt / 
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dr 



sph\ 



(9) 



At a finite temperature T the solution with period 1/T has an information on the decay 
rate[4, 13], which is a result of the saddle point approximation. 
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The possible type of transition of decay rate from thermal activity to quantum tunneling 
regimes is thoroughly discussed in Ref. [14] within quantum mechanical models. By using a 
relation 

f = * do, 

where r and E are the period of solution and energy, respectively, Ref. [14] has shown 
that the Euclidean action for the finite temperature solution meets S 3 /T-curve smoothly 
at T sph = l/r sph . The general features of action- vs-temperature for the first- and second- 
order transitions are shown at Fig.l. Although it is not general criterion for the sharp 
transition, one can see from Fig. 1(b) that S sp h is larger than S4 for the comparatively strong 
first-order transition. 

In spite of quantum mechanical ground of Ref. [14] its generalization to field theories is 
straightforward since Eq.(lO) holds even in field theories as a saddle point equation[15]. So, 
we can use a condition S sp h > S4 for a strong first-order transition. Since S3 and £4 have 
been obtained in Ref. [9] numerically 

fi4-7T 

s 3 = (2 - sy'Vsw + p 2 5 2 + (3 3 5 3 ), 

oi 

4tt 2 

S A = —(2-5)- 3 ( ai 5 + a 2 5 2 + a 3 5 3 ), (11) 
o 

where f3 1 = 8.2938, f3 2 = -5.5330, (3 3 = 0.8180, a x = 13.832, a 2 = -10.819, and a 3 = 
2.0765, the remaining one we have to calculate is r sp h. 

To obtain r sp h let us consider the small thermal fluctuation around the sphaleron: 

0(r, f) = <f) sp h(r) + ri(r, r). (12) 

Inserting it to the field equation, one can show directly that r)(T,r) obeys 

k = hr j + G 2 [ri\ + G 3 [r j ], (13) 

where 

f= 

h = -V +(30^-60^ + 5), 
G 2 [r ] ] = 3(<f )sph -l) V 2 , (14) 
Gslv]=V 3 - 



After neglecting G 2 [rj] and G 3 [rj] which are higher order terms of rj, one can show easily 

2tt 



T~sph 



h 



(15) 



where ho is negative eigenvalue of h operator. Since eigenfunction for the negative eigenvalue 
is angle-independent, we have to solve the spectra in the radial equation 

d 2 2d 



dr 2 r dr 



+ U(r) 



u n (r) = h n u n (r), 



(16) 



where U(r) = 3((f) 2 sph — 2cf) sph ) + 5. The eigenvalues for the negative mode and first positive 
mode are obtained numerically(See Fig. 2). Fig. 3 shows the numerical result of S sp h/ S± when 
Eq.(ll) is used. In this case we get a condition 5 > 5^ = 1.45 for the strong first-order 
transition. Since, however, the author of Ref.[9] obtained S3 and S4 by fitting formula with 
only three parameters, Eq.(ll) can be rough estimation. So, we calculated again S3 and £4 
numerically by using the simple shooting method, which results in 5 > 5\ = 1.42, which is 
lying within his fitting error, as shown in Fig.3. 

The opposite bound is obtained by using a criterion for the first order transition developed 
at Ref.[16]. In fact, authors of Ref. [16] obtained the condition for the sharp transition 
analytically which can be summarized as follows: 



< u I f[u ] > < 0, 



(17) 



where 



/N = - 2 *T 



rj=u 



G 2 [u ] 



(18) 



+ H, 



< a I b > is usual inner product, and u is negative mode of h. Noting that the angle 
dependent eigenfunctions do not contribute to < u \ f[u } > and using explicit forms of G 2 
and G 3 given at Eq.(14), one can reduce the condition (17) in our case to 



h + h< 0, 



(19) 



where 
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J 2 = I dk[ — + - - | - 2 j |< u fc I G 2 [«o] >l 



* 4 " ' " L " J V^o 2h + Acu 2 sph J " 1 
h = -Ji ~ h- (20) 

Here, 

J -i(^^) K|6lH>1 ' (21) 

' 1 1 

h k 2h k + 4u'i ph/ 

uo, u n (n — 1, 2, • • •), and are negative mode, discrete positive modes, and continuum states 
in the radial equation (16), respectively, and ho, h n , and hk are corresponding eigenvalues. 
Since I\ > and I2 < 0, the competition of ii and J 2 determines the type of transition. 
Numerical calculation shows that there exists only one positive discrete mode in the radial 
eigenvalue equation (16). Hence, J\ is simply 

J -fc + 5srTW) l<Ull&N>|2 - (22) 

Although the criterion for the sharp transition (19) is analytic expression, it is highly difficult 
to treat the continuum states numerically. In order to escape this difficulty we change J 2 to 
its upper bound J^ 1 
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6 26 + 40%*) 
< G 2 [u ] I G 2 [u ] > - \< u I G 2 [uo] >| 2 - |< ui I G 2 [uo] >| 



(23) 



In deriving the upper bound J^ 1 , we used the fact that the continuum eigenvalue is 5 + k 2 . 
Since J^ 1 > J 2 , h + I™ > is a sufficient condition for the second-order transition, where 
I^ 1 = — J 1 — J 2 M '■ Results of numerical calculation for Ii — J x and I\ — J\ — are given 
at Fig.4. Since h — J\ is the contributions of only negative and discrete positive modes to 
< Mo I f[uo] >, the fact that 1\ — J\ > in the full range of 5 implies the contribution of the 
continuum states are very important for the occurrence of the sharp first-order transition. 
This is new feature which does not arise at (1+1) -dimensional case, in which the contribution 
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of the continuum states is negligible [16]. Hence, we get 5 < 5 2 = 1.0245 for the second-order 
transition. Combining Fig. 3 and Fig.4 one can conclude 1.0245 < 5* < 1.42 or in terms of 
/ 0.584 < /* < 0.703. Fig.5 shows bubble nucleation at 5 = 1.18 which is calculated by 
the multigrid method, which shows the existence of the wiggly solution in the intermediate 
range of temperature. 

In this letter, we derived the lower and upper bounds of /* or 5* by applying the two 
independent criteria S sp h > S4 and I1 + I2 < to the same system simultaneously. Although 
one can obtain more accurate value of /* by solving equation of motion numerically in 
the simple model, this might be impossible for more complicated models like electroweak 
theory and real cosmological models whose potential are generally dependent on temperature 
explicitly. In these cases our method presented in this letter would be useful to determine 
the type of transition and the critical values of some parameters involved in the models. 
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FIGURES 

FIG. 1. (a) Second-order transition from the thermal to the quantum regimes, 
transition from the thermal to the quantum regimes. 



(b) First-order 



FIG. 2. Negative(Zio) and discrete positive (^i) eigenvalues of h operator. 

FIG. 3. Plot of S S ph/Si with respect to 5. The solid line is calculated by using Eq.(ll) and the 
data points are obtained by shooting method. In the range of 5 corresponding to S sp h/S4 > 1 the 
transition of the decay rate becomes comparatively strong first-order. 

FIG. 4. Plots of h - Ji and h - J ± - J 2 M with respect to 5. The fact that h - J ± > in 
the full range of 5 means that the continuum states of operator h play an important role for the 
determination of the type of transition. From the plot I± — J± — J^ 1 we can conjecture that the 
transition becomes smooth second-order when 5 < 62- 

FIG. 5. Bubble nucleation at 5 = 1.18. (a) Periodic bubble solution at T = 0.763T sp/l . (b) 
Periodic wiggly solution at T = 0.9947^. (c)0(3) symmetric sphaleron solution at high temper- 
atures. 
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Fig.1 (a) 



sph 



T 




Fig.1 (b) 



sph 



T 



2 




Fig.2 
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